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Abstract. We present numerical simulations, obtained 
with the Versatile Advection Code, of the oscillations of 
an inverse polarity prominence. The internal prominence 
equilibrium, the surrounding corona and the inert photo- 
sphere are well represented. Gravity and thermodynamics 
are not taken into account, but it is argued that these are 
not crucial. The oscillations can be understood in terms of 
a solid body moving through a plasma. The mass of this 
solid body is determined by the magnetic field topology, 
not by the prominence mass proper. The model also al- 
lows us to study the effect of the ambient coronal plasma 
on the motion of the prominence body. Horizontal oscilla- 
tions arc damped through the emission of slow waves while 
vertical oscillations are damped through the emission of 
fast waves. 

Key words: MHD - Sumcorona - Summagnetic fields - 
Sumoscillations - Sumphotosphere - Sun: prominences 



1. Introduction 

Solar prominence oscillations have been the subject of 
both observational and theoretical papers for the past 35 
years. One of the first studies (Ramsey & Smith 1966) con- 
cerned observations of global oscillations of disk filaments 
with periods of 6 m to 40 m , that were interpreted by Hyder 
(1966) as predominantly vertical motions and by Kleczek 
& Kuperus (1969) as predominantly horizontal motions. 

Tsubaki (1988) published a review on oscillation stud- 
ies of limb prominences. Most of these oscillations per- 
tain to Doppler shifts in the spectra of part of a promi- 
nence. The associated mass flows are essentially parallel to 
the photosphere, both longitudinal and transverse to the 
prominence main axis. The observed periods range from 
160 s to 82 m , with velocity amplitudes in the range of 0.2-3 
km/s. 

Zhang et al. (1991), Zhang & Engvold (1991) and 
Thompson & Schmieder (1991) studied disk filaments and 



found locally periods of 2.5 m — 22 m , with velocity ampli- 
tudes of 0.5 — 1.25 km/s. As in the observations by Ramsey 
& Smith these oscillations may have both horizontal and 
vertical components. 

Since Tsubaki's review more limb studies have been 
performed by Mashnich & Bashkirtsev (1990), Suematsu 
et al. (1990), Bashkirtsev & Mashnich (1993), Mashnich 
et al. (1993), Balthasar et al. (1993), Balthasar & Wiehr 
(1994), Park et al. (1995), Sutterlin et al. (1997) and 
Molowny-Horas et al. (1997). These observations all con- 
firm and extend previous results. There is now enough ev- 
idence to suggest the following tentative classification for 
prominence oscillations (see also Bashkirtsev & Mashnich 
1993). 

— very short periods: P w 30 s (Balthasar et al. 1993), 
perhaps due to fast waves propagating along flux tubes 
(Roberts et al. 1984). 

— short periods: P w 3 m — 10 m , at least some of which are 
related to photospheric or chromospheric forcing with 
periods of 3 m and 5 m (Balthasar et al. 1986, Zhang et 
al. 1991). 

— intermediate periods: P 10 m — 40 m , which are prob- 
ably genuine eigenmodes of the prominence (Ramsey 
& Smith 1966, Balthasar et al. 1988, Bashkirtsev & 
Mashnich 1993). Many prominences show nearly the 
same oscillation period each time they are perturbed. 

— long periods: (P w 40'™ — 114 m ), which may be related 
to chromospheric forcing (Balthasar et al. 
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We point out that the longest data set is about 7 hours 
long and the best temporary resolution is a couple of sec- 
onds. 

In particular the observed 10 m — 40 m eigenmodes are 
interesting as they are damped and thus loose energy, per- 
haps due to some interaction with the ambient corona 
(Ramsey & Smith 1966, Kleczek & Kuperus 1969). Typ- 
ically, the quality factor Q = 7rTd am p/-P < 6, with Td am p 
the damping time of the oscillation. This implies that 
three to four oscillations can be observed after the im- 
pulsive perturbation of the prominence. Understanding 
the damping mechanisms will give more insight in promi- 
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nence dynamics and may yield an extra diagnostic tool for 
prominence and ambient coronal plasma parameters. 

Prominence oscillations have been theoretically mod- 
elled by many authors. Some model the prominence as 
a harmonic oscillator with an ad-hoc damping term: Hy- 
der (1966) used viscous effects, Kleczek & Kuperus (1969) 
used emission of sound waves, van den Oord & Kupe- 
rus (1992), Schutgens (1997a) and van den Oord et al. 
(1998) used emission of Alfven waves. Other authors con- 
struct a simple MHD equilibrium and study oscillations 
thereof using the linear adiabatic MHD equations (Oliver 
et al. 1992, 1993, Oliver & Ballester 1995, 1996, Joarder 
& Roberts 1992ab, 1993, Joarder et al. 1997). Generally 
the latter approach yields marginally stable oscillations 
(no damping), although Joarder & Roberts (1992a) and 
Joarder et al. (1997) claim that leaky waves are possible 
solutions to their equations. Apparently, little effort has 
gone into studying the damping mechanisms themselves. 

Schutgens (1997ab) and van den Oord et al. (1998) re- 
cently studied the global equilibrium of prominences treat- 
ing the evolution of the magnetic field in a self-consistent 
way. The equation of motion for the prominence (ap- 
proximated as a line current) was solved simultaneously 
with the Maxwell equations for the electro-magnetic fields. 
Their results show that the Alfven travel time between 
prominence and photosphere r is an important time scale 
of the system. In particular, van den Oord et al. found 
that it is only possible to obtain stable prominence equi- 
libria by taking damping mechanisms into account. Hence, 
damping is not just necessary to describe prominence os- 
cillations quantitatively correct, but it is an essential in- 
gredient of a prominence equilibrium. 

In this paper, we study prominence oscillations, and in 
particular the effect of the ambient coronal plasma on the 
prominence motion. We use the Versatile Advection Code 
(VAC) to solve numerically the isothermal MHD equa- 
tions in two dimensions. VAC has been developed by G. 
Toth (1996, 1997) and is capable of solving a variety of hy- 
drodynamical and magneto-hydrodynamical problems in 
one, two and three dimensions using a host of numerical 
methods. In Sect. 2 a simple analytical model for promi- 
nence equilibrium and dynamics is discussed that will be 
used for comparison with our numerical results. In section 
3 we describe the isothermal equations that were solved 
numerically, the methods used, the grid structure and the 
initial conditions. In Sect. 4 the simulations are described. 
Using the analytical model mentioned before, these simu- 
lations are interpreted in terms of the physical processes 
involved. A summary and the conclusions can be found 
in Sect. 5. All dimensional variables in this paper arc in 
rational MKSA units, unless stated differently. 

2. The line current approximation 

We briefly recapitulate the line current approximation for 
prominence equilibrium and dynamics. This approxima- 



tion will serve as a guideline for discussing our numerical 
results. The prominence is approximated by a straight, 
infinitely thin and long, line current Iq > at a height 
yo parallel to the photosphere. Along the prominence we 
assume invariance. The effect of the massive photosphere 
on the prominence magnetic field is modelled through a 
mirror current —I at depth — y below the surface of the 
Sun (Kuperus & Raadu 1974, van Tend & Kuperus 1978, 
Kaastra 1985, Schutgens 1997a). The momentum equa- 
tions governing the global prominence dynamics are (ig- 
noring gravity) 



ax = - I [B, 
ay = I [B\ 



v 

mir 



(x,y) + B% OI (x,y)] - v x x, 
(x,y) + B% OT (x,y)] - v y y. 



(1) 



where a is the longitudinal mass density of the oscillating 
structure, B cor is the coronal arcade field in which the 
prominence is located and £? m i r is the field due to the 
mirror current. 

The interaction between the moving filament and the 
ambient coronal plasma gives rise to viscous effects (Hyder 
1966) and emission of magneto-acoustic waves (Kleczek & 
Kuperus 1967, van den Oord & Kuperus 1992) that act 
as damping mechanisms. These are heuristically modelled 
through v x and v y , damping constants that are in real- 
ity determined by the flow field around the prominence. 
Since viscosity of the coronal plasma is negligible, we con- 
centrate on the emission of magneto-acoustic waves. An 
approximation for the damping constants can be found 
by considering linear motion with constant velocity v of 
a solid body through a homogeneous plasma. If the cross 
section of the body perpendicular to its motion is A, the 
body transfers 2p COI cAv momentum per unit time onto the 
plasma (Landau & Lifschitz 1989, p. 256). Here c is the 
characteristic wave speed of the plasma with density p COI . 
A factor 2 is added since both the front- and backside of 
the object transfer momentum. Hence, the damping con- 
stants have the form 



v = 2p COI cA. 



(2) 



The coronal arcade is generated by a magnetic line dipolc 
M d > 0, a depth H d below the photosphere 

p M d x 2 -(y + H d ) 2 



tt {x 2 + (y + H d Y) 
2p M d x(y + H d ) 



2 - 



(3) 



7T ( X 2 + (y + H d )*y 

The mirror current's field at the location of the filament, 
in the quasi-stationary approach, is given by 

Mo la 
4-7T y ' 

Bl ir (x,y) = 0. (4) 



B^r(x,y) 



The fact that B^ ir = is a direct consequence of the mir- 
ror current mirroring the motion of the filament (to con- 
serve the photospheric flux) and the quasi-stationary field 
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assumption. When this assumption is dropped and the 
magnetic fields evolve dynamically according to Maxwell's 
laws, the expressions for B m i r become far more compli- 
cated and in particular £^ ir ^ (Schutgens 1997a, van 
den Oord et al. 1998). 

Assuming quasi-stationary field evolution (v& ~^ oo), 
prominences are in stable equilibrium provided they are 
on the symmetry axis of the arcade (y = 0) and at a height 
?/o < Ha - Furthermore, the current should attain the value 



h = 



42/oM d 
(2/o + H d f 



(5) 



Note that this prominence equilibrium has an inverse po- 
larity topology (see also Fig. 5). In fact, it corresponds to 
a Kuperus-Raadu prominence (see van Tend & Kuperus 
1978). 



where fl x and Vl y are the quasi-stationary frequencies of 
the system without damping. 

If one drops the quasi-stationary assumption, the same 
equilibrium is found, but its stability then also depends 
on the value of the coronal Alfven speed . In general, the 
solutions to the equations, which are again decoupled, are 
damped or growing (!) harmonic oscillations (Schutgens 
1997ab, van den Oord et al. 1998). 



3. Numerical methods 

We solve the time-dependent ideal isothermal MHD equa- 
tions in two dimensions using the Versatile Advection 
Code (VAC) developed by one of us (G. Toth). Descrip- 
tions of this numerical code can be found in Toth (1996, 
1997). The photosphere coincides with the y — plane. In 
conservative form the equations for density p, mass flux 
pv and magnetic field B are (p is the thermal pressure) 
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0, 

o, 
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(10) 



Note that we ignore gravity (see Sect. 5 for a discussion). 
These equations must be solved together with an equation 
of state p = c 2 p and the condition that VB = 0. Here c s is 
the sound speed, a free parameter of the system, which is 
constant throughout the numerical domain. The magnetic 
field unit is chosen such that the current density satisfies 
J = V x B (i.e. = I), ah other units are rational 
MKSA. 



Fig. 1. Location of cell centers of the grid. The photosphere 
is at 1/ = 0, and the prominence at (x,y) = (0, 3 x 10 7 ) m. 
Its radius is 1.5 x 10 7 m. The complete grid is 3 x 10 9 m by 
1.5 x 10 9 m. 



If one linearizes around this equilibrium, the equations 
of vertical and horizontal motion decouple (due to the 
symmetry of the coronal field) and both have the form of 
the familiar damped harmonic oscillator. The solutions are 
oscillations, characterized by frequencies u> and damping 
rates 5: 
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Fig. 2. Grid spacing as a function of position in the grid. The 
dashed line represents the resolution in the a;-direction, the 
dotted line in the indirection. The photosphere is at y = 0, 
and the prominence at (x,y) = (0, 3 x 10 7 ) m. 



The equations are discretized on the same grid and 
solved using a FCT (Flux Corrected Transport) scheme. 
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Fig. 3. The volume averaged horizontal and vertical momentum: pv x , y — j v pv x>y dV/V, V = 3 x 10 9 m x 1.5 x 10 9 m. Note 
the temporary increase at t = 10 000 as the artificial damping factor a is changed. 



Since there are no discontinuities in the solution, FCT per- 
forms well. Similar results can be obtained using a TVD 
(Total Variation Diminishing) Lax-Friedrichs scheme but, 
for the problem at hand, FCT is more efficient. For a 
comparison of different methods see Toth & Odstrcil 
(f996). A projection scheme (Brackbill & Barnes 1980) 
is used to keep the magnetic divergence small. Typically 
|V • B\ < J(T 3 -B/L, where B and L are characteristic 
strength and length scale of the magnetic field. 

We use a Cartesian grid of 190 x 155 cells (3 x 10° 
m by 1.5 x 10 9 m) that is strongly distorted, with the 
highest resolution at the location of the prominence and 
a much smaller resolution near the coronal edges of the 
computational domain (see Fig. f ). In a region surround- 
ing the actual prominence the grid spacing is constant. 
Typically, the grid spacing increases with 10% from cell 
to cell outside this region. As a consequence the difference 
in resolution at the prominence and at the far coronal 
edges can amount to a factor fOO (see Fig. 2). Near the 
boundaries the grid spacing is again kept constant. 

The boundary conditions are implemented using two 
layers of ghost cells around the physical part of the grid. 
The actual boundary is located between the inner ghost 
cells and the outermost cells of the physical grid. The solu- 
tion on the physical part of the grid can be advanced using 
fluxes calculated from the ghost and physical cells around 
the boundary. The prescription for the ghost cells depends 
on the specific boundary condition and may be constant 
or depend on the solution in the adjacent physical cells. 

The photosphere is much denser than the corona 
and is therefore strongly reflecting: there should be no 
mass flux or energy flux across it. We therefore choose 
the plasma density symmetric around the photospheric 
boundary (p ] 



ghost 



P 



physical \ 



and vertical momentum anti- 



pt>P hyslcal ). The other momentum 



symmetric (pv& host 
component pv x is chosen anti-symmetric as well, since 
we assume no flows along the photosphere ('no slip' con- 



dition). The magnetic field is tied to the dense photo- 
sphere and the photospheric flux (B y ) is conserved. Mag- 
netic waves should be reflected. The initial field solution 
is stored in computer memory and subtracted from the 
advanced solution. The 'linearized' field solution thus ob- 
tained is chosen to be symmetric or anti-symmetric in the 
photospheric boundary: S| host = 5f host <° + (Stoical _ 

^physical, 0^ ^ghost _ ^ghost.O ^^physical ^physical, 0^ 

In this way, one obtains a good reflection of waves for small 
perturbations. 

The other (coronal) boundaries do not coincide with 
a physical boundary and should be as open as possible. 
Since any choice of boundary condition will always gener- 
ate some reflection, which we want to avoid, we decided to 
place these boundaries at large distances from the promi- 
nence so that twice the wave crossing time (prominence- 
coronal boundary) is longer than the simulation time. 
Also, the coarsening of the grid damps the outward mov- 
ing waves and the reflection is minimized. At the location 
of the coronal boundaries we prescribe fixed B ghost , but 
copy p and pv from the physical part of the grid to the 
ghost cells. 

The initial configuration is a superposition of three 
different analytical MHD equilibria. The global coronal 
structure is a potential arcade, given by Eq. (3). Since we 
ignore gravity, the plasma density in this arcade is con- 
stant. To this equilibrium we add the fields and plasma 
densities of two current carrying flux tubes, one above, 
the other below the photosphere. Both flux tubes are lo- 
cated on the polarity inversion line of the arcade. The flux 
tube at a height yo above the photosphere represents the 
prominence. Its total axial current is Iq. The flux tube 
—yo below the photosphere represents the mirror current 
and has a total current — Jo- 

The flux tube equilibrium is derived starting from a 
current profile (see also Forbes 1990) 
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'jo r<r a -Ar , 
3z(r) = I Jo sin 2 (f ^) r - Ar < r < r , (11) 

r > r . 

The radius of the tube ro is larger than Ar , the size of 
the region in which the current density drops off to zero. 
The total axial current in the flux tube is 



rr 

I = 2n r'j z (r')dr' 
Jo 



(12) 



and the azimuthal field component follows from Stokes' 
theorem 



2frr 



r < r , 
r > r . 



(13) 



The current is confined to r < r . Beyond this range, the 
field is potential. The pinching field B$ has to be balanced 
by gas pressure 



p(r) = - [ B^r'^r) dr, 

J rn 



(14) 



where it is assumed that the flux tube proper carries no 
mass outside r > r . The longitudinal mass density of the 
flux tube equals 



/■ro 

a a = 2n r'p(r')/c 2 s dr' , 
Jo 



(15) 



while the total longitudinal mass density of the promi- 
nence is (due to the superposition of equilibria) 



a = (T + np COT r Q . 



(16) 



In a global sense, equilibrium is obtained when Eq. (5) is 
satisfied. However, especially near the outer edge of the 
prominence no force balance exists. It is necessary to let 
the initial configuration relax to a numerical equilibrium, 
before studying oscillations. As it turns out, force imbal- 
ance is mostly due to the gradients in density and field 
within the flux tube (i.e. numerical discretization errors) 
and a stable equilibrium is readily obtained. To prevent 
the flows at the outer edge of the prominence from be- 
coming too large during the relaxation phase, and cause 
an instability, an artificial damping term dtpv = —apv is 
used. This term is switched off during subsequent oscilla- 
tion studies. 

Once a numerical equilibrium has been found, we per- 
turb it and study the resulting oscillations. In all cases, 
the perturbation is caused by instantaneously adding mo- 
mentum to the prominence. 

To check the validity of our results, we performed os- 
cillation simulations for the same physical parameters, but 
different numerical parameters. In particular we changed 
the grid resolution (Air = 5 x 10 5 or 10 6 m inside the 
prominence), the amount of stretching of the grid (0%, 
5% and 10%), the duration of the relaxation (T = 9 000 



lxio 



S -1X10 



< 



-2x10 



-3x10 




5.0xl0 3 l.OxlO 4 1.5xl0 4 2.0xl0 4 



time (s) 



Fig. 4. The change in height of the prominence during relax- 
ation. Initially the prominence was at yo = 3 x 10 7 m. During 
relaxation the prominence moves downward a mere 0.1 % of 
its initial height. 



or T = 20 000 s) and the constraint on V • B\ (< 
10~ 4 , IQ~ 3 B/L) and found that the results are similar. 
Using 0% stretching implies that the grid is rather small in 
its physical size, due to computational limitations. Reflec- 
tion at the coronal boundaries will then limit the useful- 
ness of the simulations to the first 4000 s for p cor = 10~ 12 
kg/m 3 . 





Fig. 5. The magnetic field topology after relaxation. A strong, 
localized current runs near (x,y) = (0, 3 x 10 7 ) m, at the lo- 
cation of the prominence. 



In addition, the effect of the photospheric boundary 
condition was studied. A perfectly reflecting photosphere 
does not allow a net momentum or energy flux. Also, the 
magnetic flux distribution is constant. The implementa- 
tion of the photospheric boundary conditions automati- 
cally ensures zero momentum flux and constant magnetic 
flux. However, since the magnetic field is not symmetric 
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Fig. 6. The plasma-/? (a) and Alfven speed (b) in and near the prominence. In the left graph, gray denotes f3 > 1. The solid 
line is the contour line f3 = 0.1. In the right graph, darkest gray denotes v\ < c s and lightest gray denotes va > 500 km/s. 
Maximum Alfven speed is 873 km/s. The dotted circle represents the boundary of the prominence, the dotted rectangles the 
coronal surface over which the Poynting flux was integrated. The graphs are valid for p COI — 10~ 12 kg/m 3 . For other coronal 
densities the scalings /3 and va oc l/^/p cor can be used. 



in y = 0, a small Poynting flux is present. Consider a 
box-shaped surface around the prominence. The coronal 
'walls' of this surface each have a shortest distance to the 
prominence center of 2.5 x 10 7 m (see Fig. 6) and reach 
down to the photosphere, along which we close the box. 
The total time-averaged Poynting flux through the pho- 
tosphere is typically 1CT 3 — 10~ 2 times smaller than the 
coronal flux, for vertical oscillations. For horizontal oscil- 
lations, the coronal flux is small in essence (due to the 
anti-symmetry) and provides no reasonable yardstick (of 
course one could still use the flux through a single 'wall', 
but this is bound to lead to similar results as for vertical 
oscillations). Furthermore, we point out that the equilib- 
rium height of the prominence during oscillation studies 
does not change by more than 0.06%, and often even less. 
We therefore feel that the photospheric boundary is suffi- 
ciently well represented in our numerical boundary condi- 
tions. 

4. Oscillation studies 

Before we study oscillations in prominences, we first have 
to compute a stable, numerical prominence equilibrium. 

For the coronal arcade, we choose the following param- 
eter values: = 4 x 10 s m, Md = 10 20 Am and plasma- 
density /9 cor = 10~ 12 kg/m 3 . The two flux tubes are given 
by: r = 1.5 x 10 7 m, Ar = 10 7 m and j = 2x 10~ 4 
A/m 2 . Hence the longitudinal mass density of the flux 
tube proper is cr — 1.28 x 10 4 kg/m, and the density of 
the prominence is a = 1.35 x 10 4 kg/m. The total current 
is Iq = 6.5 x 10 10 A. Such a prominence is expected to 
oscillate with a horizontal period of 2777 s and a vertical 
period of 1118 s. The sound speed was chosen c s = 128.5 
km/s, typical of the corona at a temperature T = 10 6 K. 



We now let the initial state as defined in the previ- 
ous section relax to a numerical equilibrium. The artifi- 
cial damping term (dtpv = —apv) was used to prevent 
numerical instabilities. We choose a = 0.1 for the first 10 4 
seconds of the relaxation, and a = 0.01 during the latter 
10 4 seconds. Figure 3 shows the relaxation of the volume 
averaged horizontal and vertical momentum. At the end 
of the relaxation the flows in the larger part of the arcade 
are typically less than 5 m/s. Along the boundary of the 
prominence body a rather irregular flow field exists with 
velocities of 175 m/s at most. The prominence still moves 
downwards at a systematic speed of some 1 m/s (see Fig. 
4). Considering both the magnitude of the velocity pertur- 
bation applied later and the total simulation time scale, 
we consider these residual flows to be unimportant. The 
agreement between simulations starting from different re- 
laxations (with a duration of either 9 000 or 20 000 s) con- 
firm this. 

The numerical equilibrium does not deviate much from 
the initial state. The resulting field topology around the 
prominence is shown in Fig. 5. The prominence itself is a 
region of high plasma- /3, surrounded by a region of mag- 
netically dominated plasma (see Fig. 6). At large distances 
from the prominence the plasma-/? becomes larger than 
unity again, as the magnetic field strength decreases while 
the coronal density is constant (due to the absence of grav- 
ity). The Alfven speed is also shown in Fig. 6. It increases 
as one approaches the prominence from the corona, once 
inside it falls off rapidly. The sound speed, of course, is 
the same everywhere (it is a free parameter of Eq. (10)). 

The obtained numerical equilibrium is used to derive 
a series of numerical equilibria by adding (or subtracting) 
a constant value from the density in each cell. Since this 
does not create any additional forces (the force due to 
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Fig. 7. The volume averaged horizontal and vertical momenta pv x , y as a function of time for all eight simulations. Curves 
with a higher initial value of the momentum correspond to larger coronal densities. All simulations were started with the same 
velocity perturbation (10 km/s). For the sake of clarity only the first 5000 s are shown for the vertical momentum. 



Table 1. Characteristic time scales of prominence dynamics for all eight cases considered. The coronal density p COI and (as a 
consequence) the prominence longitudinal mass density a are the only parameters that differ. In particular the total current 
Jo = 6.5 x 10 10 A and the initial perturbation (10 km/s) are the same for all simulations. The estimated error in the period 
is typically less than 3%, vertical periods having larger errors than horizontal periods. The estimated error in the damping 
time is typically less than 10%, vertical damping times having smaller errors than horizontal damping times. For the horizontal 
oscillations, we used only cases 1-6, for vertical oscillations only cases 1-5 to compare the numerical results to the simple line 
current model of Sect. 2. No reliable vertical time scales could be obtained for the last two cases due to strong transients and 
noise. 
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io- 13 


1.29 x 


10 4 


1.47 x 


10 4 


2 961 


26 600 


1 118 


1100 


3 


5 x 


io- 13 


1.31 x 


10 4 


1.68 x 


10 4 


3134 


16100 


1087 


860 


4 


1 X 


io- 12 


1.35 x 


10 4 


2.13 x 


10 4 


3 461 


9 400 


1053 


820 


5 


2 x 


io- 12 


1.42 x 


10 4 


3.07 x 


10 4 


4160 


5 600 


1080 


990 


6 


4 x 


io- 12 


1.56 x 


10 4 


5.04 x 


10 4 


5 600 


3 800 


1350 


1020 


7 


8 x 


io- 12 


1.84 x 


10 4 


9.20 x 


10 4 


7 740 


8100 


?? 


?? 


8 


1.6 x 


io- 11 


2.41 x 


10 4 


17.95 x 


10 4 


10 300 


22 200 


?? 


?? 



gas pressure is — Vp) a new equilibrium is found. In the 
resulting equilibria we disturb the inner part of the promi- 
nence (all plasma within 1.25 x 10 7 m from the center of 
mass) with a velocity perturbation of 10 km/s, at an angle 
of 45° to the photosphere. The evolution of the system is 
studied for 10 4 s, in some cases even for 2 x 10 4 s. 

In Fig. 7 the evolution of the volume averaged mo- 
mentum is shown for all eight cases considered. Initially 
momentum is concentrated in the prominence, but it is re- 
distributed throughout the corona in the subsequent evo- 
lution. A global oscillation is apparent, whose properties 
depend strongly on coronal density. 

For one particular case, the flow field is shown in Figs. 
8 and 9. Although only the velocity is shown, the promi- 
nence stands out clearly in most graphs. It seems to move 
through the coronal plasma as a rigid body. The largest 



velocities are usually found outside the prominence. The 
coronal plasma 'washes around the filament'. 

In particular, we studied the motion of the prominence. 
This was done by computing every so many time steps 
the center of longitudinal mass and current density. To 
rule out the contribution of the coronal part of the grid, 
only densities above a certain threshold where used. As 
the contrast between typical coronal and prominence cur- 
rent densities is larger than the contrast between typical 
coronal and prominence mass densities, the first provide 
a better estimate of the location of the prominence. Nev- 
ertheless, the results are always similar. 

In essence, the motion of the filament can be described 
by two decoupled damped harmonic oscillators. The hor- 
izontal resp. vertical displacement of the center of longi- 
tudinal mass or current density of the prominence may 
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Fig. 8. Flow field of the plasma in the prominence and ambient corona, during the first f 600 s. Shown are the numerical results 
for p C or = 10~ 12 kg/m 3 , see Table 1. Maximal velocity is 11.7 km/s. 



be fitted to Ae~*/ Tdamp sm2nt/P. The resulting periods 
and damping times for the horizontal and vertical oscil- 
lations are listed in Table 1. Sometimes strong transient 
effects at the start of the simulation and noise at later 
times (when the velocities are smaller) cause deviations 
from a simple damped harmonic oscillator. The horizon- 
tal oscillations yield, in general, better fits. Typically the 
low density simulations lead to better fits than the high 
density simulations. For horizontal motions we used the 
cases 1-6, for vertical motions we used the cases 1-5 (see 
Table 1). No reliable vertical time scales could be obtained 
for the last two cases due to strong transients and noise. 

Comparing the results with simulations where a purely 
vertical or horizontal perturbation was applied shows that 
the horizontal and vertical motions of the filament are ac- 
tually decoupled. This is understandable as the filament 
is located on the symmetry line of the coronal field (see 
also van den Oord et al. 1998). At the same time this sug- 
gests that non-linear effects are not important. Further 
evidence for the linearity of the numerical results is found 
in the good fits of the prominence oscillation curves to 
a damped harmonic oscillator. However, simulations with 
smaller or larger perturbations (3 km/s or 20 km/s in- 
stead of 10 km/s) yield different damping times for the 
horizontal motions (see Fig. 10). For all three simulations 



the horizontal periods differ by only 0.1%, while the hori- 
zontal damping times differ by 30 — 60%! Apparently, the 
horizontal damping mechanism depends non-linearly on 
the flow speeds. We have no explanation for this result, 
but believe it to be genuine. Extensive tests with higher 
spatial and temporal resolution argue strongly against the 
possibility of a numerical effect. It is a strong indication 
that the damping mechanism for horizontal and vertical 
motion differ, at least in their dependencies on the flow 
field. 

The horizontal periods show a strong dependence on 
the coronal density, while the vertical periods are almost 
constant (Fig. 11). The quality factors (Fig. 12) of the hor- 
izontal oscillation are typically larger than four, and hence 
damping contributes little to the oscillation frequency (see 
Eq. (6)). The dotted lines (in Fig. 11) represent the un- 
damped quasi-stationary period for a prominence with a 
longitudinal mass density of a = ctq + ttvq p cor . Clearly this 
is a bad fit to the numerical results for horizontal oscilla- 
tions. Let us assume that, due to the magnetic structure, 
the actual oscillating body has a radius r e s > r$. The 
actual oscillating body now has an effective longitudinal 
mass density a c g = cr Q + 7rr^ ff p cor . From Eq. (6) we obtain 

Pcor 



2/o 



^ (2/o + H d f 



P z 



<?0 



(17) 
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Fig. 9. Flow field of the plasma in the prominence and ambient corona, during the latter part of the evolution. Shown are the 
numerical results for p COI — 1CP 12 kg/m 3 , see Table 1. Maximal velocity is 8.5 km/s. 



the left-hand-side of which can be fitted to a func- 
tion of the form Cp2 OI - We find that r e ff = 5.0 x 
10 7 (jO C or/10~ 12 ) ' 04 m, weakly dependent on the coronal 
density. The typical radius r e s = 5.0 x 10 7 m agrees well 
with the location of the magnetic surface across which the 
connectivity of the field lines changes from closed field 
lines in the corona to field lines anchored in the pho- 
tosphere. An even better approximation to this surface 
is obtained by taking into account the effect of induced 
mass as described by Landau & Lifschitz (1989, p. 29; 
see also Lamb 1945). For an incompressible, potential hy- 
drodynamical flow the effective mass is the sum of the 
mass of the actual oscillating body (er + nrlp COI ) plus 
the mass of the fluid displaced by the body. In that case 
c e ff = °"o + 27rr 2 ff/ o cor (note the factor 2!) and we find 
r cff = 3.5 x 10 7 ( Pcor /10- 12 ) - 04 m. 



In Fig. 13, we have plotted the ratio of vertical damp- 
ing time to horizontal damping time. From Eqs. (2), (7) 
and (9) this ratio is proportional to c x /c y . Here c x is 
the typical wave speed for perturbations travelling par- 
allel to the photosphere and c y is the typical wave speed 
for perturbations travelling perpendicular to the photo- 
sphere. For simplicity, we assume that the cross sections 
for both directions have the same p COI dependence, but 
given the azimuthal invariance of the prominence flux tube 
this is probably a fair approximation. Now c x and c y are 
equal to either the slow cusp speed op, the Alfven speed 
ca or the fast speed c t . In a low plasma-/? environment, 
the cusp speed equals the sound speed c s , while the fast 
speed equals the Alfven speed ca- Hence c x /c y cx p^"' 5 , 1 
or p2or depending on whether slow, Alfven or fast wave 
emission prevails in a certain direction. From the data, we 
find c x /c y oc p^' 53 which strongly suggests that slow wave 
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Fig. 10. The center of the current density above the threshold 
value. Dashed line: 20 km/s perturbation; solid line: 10 km/s 
perturbation; dotted line: 3 km/s perturbation. The first curve 
was multiplied by 0.5, the last curve was multiplied by 3.333 
to facilitate comparison. The top graph shows the horizontal 
displacement, the bottom graph the vertical displacement. The 
coronal density equals p cor = 2.5 x 10~ 13 kg/m 3 . 



Fig. 12. Quality factors (Q — 7rTdam P /-P) for the horizontal 
(squares) and vertical (diamonds) oscillations. 
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Fig. 11. Periods for the horizontal (squares) and vertical (di- 
amonds) oscillations. The dotted lines are the quasi-stationary 
periods for a prominence loaded with mass a = <r + Tvrop COI . 
The top dashed line is the model for horizontal oscillations as- 
suming that the actual oscillating body is larger than just the 
prominence proper. The bottom dashed line is the model for 
vertical oscillations using the effective longitudinal mass den- 
sity <T e ff as determined from the horizontal oscillations. 



emission damps horizontal motions, while Alfvcn or fast 
wave emission damps vertical motions. Since the promi- 
nence moves vertically to the field lines in the latter case, 
fast waves are more likely. 

The conclusions regarding damping mechanisms can 
be substantiated further by comparing the relevant time 
scales for horizontal and vertical oscillations indepen- 
dently to the Eqs. (7) and (9) (see Fig. 14). For the mass, 
we use the effective longitudinal mass density a e g as ob- 



density (kg/m 3 ) 



Fig. 13. Ratio of the damping times of vertical and horizontal 
oscillations. The solid line (oc p^r' 53 ) is a fit to the data points 
(plus signs). The dotted lines are fits to c x /c y oc pcor 5 , 1 or p cor 
using the data point for p cor = 5x 10~ 13 kg/m 3 as a gauge. 



tained previously. From Eq. (2) and either Eq. (7) or (9) 
we find 

= Ac, (18) 

Pcor-' damp 

the left-hand-side of which can be fitted to Cp2 OI - Here 
A is the cross-section perpendicular to the direction of 
motion. For horizontal oscillations we find c x oc Pcor 6±0 ' 06 : 
and A w 2 x 10 7 m 2 . This suggests excitation of slow 
waves by horizontal motion of the prominence. For vertical 
oscillations c y oc Pc^ 54 ^' 04 which suggests excitation of 
fast waves by vertical motion of the prominence. Since the 
Alfven speed varies in space, it is not possible to derive a 
value of A from the product of Ac. 

The cross-section A for the horizontal damping mech- 
anism yields a dimension for the oscillating body smaller 
than the extended size r c ff as obtained from the horizontal 
periods. This may be due to field line curvature. Since the 
size of the oscillating 'solid' body is determined by the field 
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Fig. 14. The effective longitudinal mass density a e ff over 
the damping time Td and coronal density p CO i for horizontal 
(squares) and vertical (diamonds) oscillations. The dotted lines 
are model fits to Eq. (18). 

topology, Lorentz forces cause the waves that carry away 
momentum. However, Lorentz forces only act perpendic- 
ular to, not parallel to the field lines. Thus only perturba- 
tions of magnetic field with a strong vertical component 
can excite horizontally travelling waves. 

5. Summary and conclusions 

We have made a numerical investigation of prominence 
oscillations, by solving the isothermal MHD equations in 
two dimensions. First we computed a prominence equilib- 
rium that is very similar to the Kuperus-Raadu (inverse 
polarity) topology. However, in our numerical equilibrium 
the prominence is not infinitely thin, but instead well re- 
solved. From this equilibrium, we derived other equilib- 
ria with different coronal plasma densities. We then per- 
turbed the system by instantaneously adding momentum 
to the prominence mass and followed the ensuing oscil- 
lations. The dependence of the characteristic time scales 
(periods and damping times) on the coronal plasma den- 
sity was analyzed in terms of a solid body moving through 
a fluid. To our knowledge this is the first attempt at nu- 
merically simulating prominence oscillations. 

In our numerical model, we ignored the effects of grav- 
ity and thermodynamics, for the sake of clarity and practi- 
cality. Also, up to now the numerical boundary conditions 
we are using do not seem to allow a stable gravitationally 
stratified corona. 

We believe, however, that gravity and thermodynamics 
do not contribute significantly to the physics of the sys- 
tem. Gravity is of small consequence for the equilibrium 
of a Kuperus-Raadu prominence as detailed by van Tend 
& Kupcrus (1978). Force balance is due to two Lorentz 
forces: one due to the coronal magnetic arcade, the other 
due to the photosphcric flux conservation. The gravita- 
tional force can be ignored when describing global equilib- 



rium. Furthermore, the scale height of a corona of T = 10 6 
K is 3 x 10 8 m, which is larger than the typical vertical 
size of prominences. The inclusion of gravity would change 
the appearance of the prominence into a slab, with promi- 
nennce matter accumulating in the pre-existing (!) dips in 
the magnetic field lines. The magnetic field configuration 
would hardly change. 

Likewise, the absence of a thermally structured corona 
and prominence seems of only minor influence. The most 
important effect would be that a cool prominence will be 
heavier than a prominence of coronal temperatures. In 
our model this is balanced by the absence of a longitu- 
dinal field. As a consequence, the pinching effect of the 
azimuthal prominence field has to be balanced by gas pres- 
sure only (in real prominences the longitudinal field pres- 
sure contributes significantly) . The total longitudinal mass 
density of our prominence (~ 1.3 x 10 4 kg/m 3 ) agrees well 
with that of a slab of height 5 x 10 7 m, width 6 x 10 6 m and 
density 5 x 1CT 11 kg/m 3 . Oliver & Ballester (1996) stud- 
ied the influence of the prominence-corona transition re- 
gion (characterized by strong temperature gradients) and 
found that it mainly influences the prominence internal 
oscillations, but not the global oscillations. 

We point out that the inclusion of both gravity and 
thermodynamics would give the prominence proper the 
appearance of a cool slab, suspended in the dips of the 
field lines belonging to the prominence current. Presum- 
ably this will not change the global oscillation discussed 
in this paper, since that mode is determined by the overall 
field structure. 

The results indicate that, for typical coronal densities 
(pcor = 10 -13 — 1CP 12 kg/m 3 ), the prominence structure 
can indeed be viewed as a solid body moving through a 
fluid. However, the mass of this solid body is determined 
by the magnetic topology, not the prominence proper. In 
particular, the mass of the body seems to be determined 
by coronal field lines that enclose the prominence proper. 
In a low plasma-/? environment this is to be expected. As 
a consequence, the total mass of the oscillating solid body 
is larger than the mass of the prominence proper. 

Due to the symmetry of the system, horizontal and 
vertical prominence oscillations decouple. These oscilla- 
tions can each be interpreted as the motion of a damped 
harmonic oscillator. The horizontal periods and the hor- 
izontal and vertical damping times can be explained by 
assuming that the actual oscillating structure is larger 
than the prominence proper, due to the magnetic field 
topology. Only the vertical periods do not agree with this 
model. They are nearly constant for different values of 
the coronal density and are best modelled by the oscilla- 
tion of the prominence proper in a corona with vanishing 
plasma density (p COI J. 0). For realistic coronal densities, 
the horizontal periods change with 3 — 15% at most when 
taking the actual oscillating structure into account. How- 
ever, for lower prominence longitudinal mass densities a 
in a stronger coronal background field (and hence larger 
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prominence currents Io), the effect will be much more pro- 
nounced. 

Vertical oscillations lead to the emission of fast waves 
that carry momentum away from the prominence and 
damp the oscillation. Horizontal oscillations, on the other 
hand, lead to the emission of slow waves. These will damp 
the horizontal oscillation, but less effectively than fast 
waves (Q x > Q y ). 

The difference in wave emission between horizontal 
and vertical oscillations can be understood in terms of 
the coronal arcade in which the prominence is embed- 
ded. Due to the large scale height, the arcade field is close 
to horizontal near the prominence. Waves that travel in 
the vertical direction (up or down) therefor travel more 
or less perpendicular to the field lines and must be fast 
waves. Waves that travel in the horizontal direction travel 
along the field lines. They could be either Alfven waves 
or magneto-acoustic waves. As regards excitation of the 
waves, it is obvious that the vertical motion of a promi- 
nence across field lines that are nearly horizontal will com- 
press both gas and magnetic field and thus set off fast 
waves. The excitation of the magneto-acoustic slow waves 
(for our analysis strongly suggests they are slow waves) 
is not as well understood. But apparently the prominence 
acts as a piston during horizontal motions along magnetic 
field lines of the arcade. 

We surmise that for smaller scale height of the arcade 
the slow waves might be replaced by fast waves. 
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